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Taste breaking in staggered fermions from random matrix theory 
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We discuss the construction of a chiral random matrix model for staggered fermions. This model includes 
0(a 2 ) corrections to the continuum limit of staggered fermions and is related to the zero momentum limit of the 
Lee-Sharpe Lagrangian for staggered fermions. The naive construction based on a specific expansion in lattice 
spacing (a) of the Dirac matrix produces the term which gives the dominant contribution to the observed taste 
splitting in the pion masses. A more careful analysis can include extra terms which are also consistent with the 
symmetries of staggered fermions. Lastly I will mention possible uses of the model including studies of topology 
and fractional powers of the fermion determinant. 



Staggered fermions are a popular way to in- 
clude quarks on the lattice. A single staggered 
fermion determinant in the continuum limit ac- 
tually describes four tastes of quarks, so named 
to distinguish them from physical flavors. At fi- 
nite lattice spacing the tastes interact and the 
continuum SU (4) taste symmetry is broken. The 
size of the taste breaking can be reduced by using 
improved actions, but for current lattice spacings 
the effects are clearly present in chiral fits pp. 

The form of the taste breaking terms for a sin- 
gle flavor (four tastes) of staggered fermions was 
worked out as corrections to the chiral Lagrangian 
by Lee and Sharpe ■ They included all terms of 
0(a 2 ) that are consistent with the symmetries of 
staggered fermions. This was later extended to 
arbitrary flavor content (most notably 2+1 fla- 
vors) by Bernard and Aubin 3 . 

Here we will outline how to include taste break- 
ing in a chiral random matrix model. This model 
should be equivalent to low energy lattice QCD 
with staggered fermions and can be used to study 
properties of the lattice theory which are sensitive 
to the smallest eigenvalues of the Dirac matrix. 
A more detailed study of the staggered matrix 
model and its applications will appear elsewhere. 



1. Random Matrix Theory 

Chiral random matrix theory has been used 
extensively to model the low energy eigenmodes 
of the Dirac operator. Its construction is based 



solely on the symmetries of the Dirac operator. 
In a chiral basis the standard Dirac operator for 
a single quark flavor takes the form 
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where W is a complex N x N matrix, In is the 
N x N identity matrix and m is the quark mass. 
Here we will only consider the case of zero topo- 
logical charge and three or more colors which cor- 
responds to the Dirac matrix given above. Due to 
universality the partition function near the chiral 
limit is not sensitive to the exact form for the dis- 
tribution of the various matrix elements. One can 
therefore make the simplest choice of an indepen- 
dent Gaussian measure for each matrix element 

H(W) = exp {-a N Tr W f W ) . (2) 

The parameter a is related to the chiral limit of 
the chiral condensate (S) and the volume (V) by 
i/a = VE/2N. For a recent review of chiral ran- 
dom matrix models see 0] and its references. 

In order to extend chiral random matrix the- 
ory to include taste breaking we need to examine 
the form of the taste breaking interactions. The 
original motivation for this work came from ex- 
amining the classical expansion of the staggered 
Dirac operator near the continuum limit given by 
Kluberg-Stern et. al. 0. Their expansion starts 

(7 M <X> I 4 ) + (I 4 <X> h) m + a (75 <g> 6,5) D 2 V . (3) 

The notation (S <g> T) has the standard meaning 
of the outer product of a 4 x 4 spin matrix S with 
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a 4 x 4 taste matrix T with £„ 6 = £„£ 5 = 7*75. 
The first part is just the standard Dirac oper- 
ator for four identical flavors of mass m. The 
remaining term is suppressed by a factor of the 
lattice spacing and breaks the SU(A) taste sym- 
metry. There is also an 0(ag) term which for our 
purposes is similar to the included taste breaking 
term and will be ignored. The remaining correc- 
tions are at least 0(a 2 ). This form of the expan- 
sion is specific to the particular basis chosen for 
the continuum Dirac fermions. In fact Luo |B] 
has shown that with an improved basis the 0(a) 
terms can be removed giving corrections only at 
0(a 2 ). This ambiguity in basis to use for the 
continuum expansion can be avoided by instead 
starting with the same four-fermion operators in 
the continuum effective Lagrangian used by Lee 
and Sharpe 0- Using a Hubbard-Stratonovitch 
transformation these terms can split into fermion 
bilincars which can be used to construct a matrix 
theory with the correct symmetries. We do not 
present the details of this construction here but 
will refer to it later in the paper. 

For simplicity here we will continue to con- 
struct the matrix model based on the expansion 
© that was given. The taste breaking term in 
the chiral basis has the form 

T = (A u <g> cr 3 ) ® &5 (4) 

where A v is an N x N Hermitian matrix and 03 a 
Pauli matrix. We again make the simplest choice 
for the integration measure over the matrices A v 
which is that of independent Gaussian matrix el- 
ements given by 

fJ-(A u ) — exp (—(3 NTy A 2 ) (5) 

with some unknown parameter p. The partition 
function is then 

Z = I (j,(W)dW Y[n(A v )dA v Bet(M) (6) 
with 

M = V ® I 4 + a T . (7) 

We now will relate this model to the staggered 
chiral Lagrangian given by Lee and Sharpe. 



2. Chiral Lagrangian 

The random matrix model © can be trans- 
formed into a nonlinear sigma model using stan- 
dard methods. In the large N limit the sigma 
model can be expanded around its saddle point 
which gives 

/ dU exp(-mU£ReTrt/ - a 2 V ) (8) 

J 17(4) 

with 

V = N 7iYl KeTr (UUUfa,) ■ (9) 

This is equivalent to the zero momentum limit of 
an effective chiral theory. Notice that the taste 
breaking is 0(a 2 ) in the chiral Lagrangian even 
though we started with a term of 0{a) in the 
Dirac operator. The form of V is identical to 
the term with coefficient C4 in the Lee-Sharpe 
Lagrangian if we set Not/ (3 — VC4. It is in- 
teresting to note that this term was previously 
found to be the dominant contributor to the mass 
splittings between pions with different taste struc- 
tures. In FigureQlwe show the difference between 
the squares of various pion masses and the square 
of the Goldstone pion (£5) mass in units of the lat- 
tice spacing. On the x-axis are the taste genera- 
tors used with the corresponding contribution to 
the splitting from the C4 term written below. The 
data are from simulations of improved "a 2 -tad" 
staggered fermions at (3 imp — 7.3 and am — 0.02 
tabulated in [7] . The results agree very well with 
the linear rise predicted from the C4 term. 

We have so far considered only the taste break- 
ing term given in 0. A more general construc- 
tion alluded to earlier gives the following addi- 
tional terms for the matrix model 

(Ai 1] <g> 00,3) -> C* 3 

(A (2) <g) <ti) ® H4, — > const. 

(A< 3 > «) cti) <g> i& -> Ci 

(A^®^)®^ C 6 . (10) 

The (To, 3 can be either <7o = I2 or 173. Also the 
03 in (0J can be replaced with <ro to produce an- 
other valid term. The arrow shows which term 
from the chiral Lagrangian the matrix model term 
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Figure 1. Pion mass splittings for different taste 
content relative to the Goldstone pion. The con- 
tribution from the C4 term of the staggered chiral 
Lagrangian is given at the bottom. 

produces in the large N limit. The corresponding 
terms in the chiral Lagrangian are 

* — 'm 

CeE TrttuvU^tf). (11) 
The two remaining terms from [2] are 

C 2 Rc[Tr([/ 2 ) - TrU 5 £/£ 5 [/)] (12) 
C 5 ^V Tr(^L/^[/t)-Tr(^ 5 [/£ 5AI [/t). 

They can also be produced from the matrix model 
using terms similar to those in (|10l) , however some 
of the terms would have to be Hermitian (note 
the terms in H10|) are all anti-Hcrmitian) which 
would make them not suitable for modeling the 
Dirac matrix. This point requires further study. 

3. Conclusions and Future Applications 

We have presented a chiral random matrix 
model for staggered fermions. In its simplest form 
© it includes just the term that is known to be 
dominant in the pion mass splittings. Additional 



terms consistent with the staggered fermion sym- 
metries can also be included, but the construction 
of the most general form of the matrix model still 
requires further study. 

This model can also be extended to include ad- 
ditional flavors and non-zero topological charge. 
Both of these extensions will be covered in a fu- 
ture publication. Since topology is related to the 
lowest eigenvalues of the Dirac matrix, this should 
prove to be a useful model for studying the topol- 
ogy of staggered fermions. 

Current lattice simulations using staggered 
fermions typically take the square root or fourth 
root of the fermion determinant to reduce the 
number of tastes from four to produce one or two 
physical flavors in the continuum limit. The ef- 
fects of the root on the low energy spectrum are 
not well understood. One could study its effects 
in the matrix theory where analytic results are 
easier to produce. If the distribution of Dirac 
eigenvalues in the matrix model agrees with lat- 
tice calculations, then the effect of the root should 
be modeled properly by the matrix theory. This 
is currently perhaps the most interesting applica- 
tion of the staggered random matrix model. 

I would like to thank Claude Bernard, Car- 
leton DeTar and Urs Heller for discussions which 
helped me understand taste breaking in staggered 
fermions better. This work was funded by by the 
US NSF. 
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